Probability current tornado loops in three-dimensional scattering 
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We consider scattering of a three-dimensional particle on 
a finite family of 8 potentials. For some parameter values 
the scattering wavenctions exhibit nodal lines in the form 
of closed loops, which may touch but do not entangle. The 
corresponding probability current forms vortical singularities 
around these lines; if the scattered particle is charged, this 
gives rise to magnetic flux loops. The conclusions extend to 
scattering on hard obstacles or smooth potentials. 



The fact that quantum systems may exhibit nontrivial 
topological effects is known for long [Q. Recent interest 
to vortices has been focused mostly at superconducting 
systems described by the Ginzburg-Landau equation . 
However, vortices have been also observed in pure quan- 
tum mechanics, specifically in numerical analysis of var- 
ious models of mesoscopic electron transport |3j where 
they may arise even without an applied magnetic field. 

The vortical behavior in quantum mechanical scatter- 
ing is closely relate d with the wavefunction phase. Writ- 
ing = \J p(r) e l< ^ r ) , we can express the probability 
current as j(f) = yo(r)V(/)(r) Q|. In a region where ex- 
ternal forces are absent, the integral of j(r) over a closed 
loop can be nonzero only if it encircles a singularity in 
which the phase 4> is ambiguous. However, solutions 
to the stationary Schrodinger equation are smooth func- 
tions, so only such singularities are zeros of tp . 

The existing results mentioned above are usually con- 
cerned with two-dimesional systems where the vortices 
are planar and centered around nodal points. The aim of 
the present letter is to show how the probability flow vor- 
tices look like in three dimensions, where nodal sets of the 
scattering wavefunctions are generically smooth curves. 
While the examples of Ref. H deal with perturbed chan- 
nels of various shapes, the presence of the boundary is in 
fact not essential for the existence of vortical solutions. 
To illustrate this we shall discuss scattering in Ft 3 . For 
the sake of simplicity we are going to analyze a simple 
example of a particle scattered on a finite family of point 
interactions, however, the conclusions extend to a much 
wider family of potentials. 

We shall demonstrate that scattering wavefunctions 
have often nodal lines, and that the latter are of the form 
of closed loops. The probability current in the vicinity of 
the line is locally cylindric (tornado-shaped). Moreover, 



if the scattered particles are charged, the corresponding 
electric current generates a magnetic field whose flux lines 
not far from the nodal lines are closed. 

Let us describe the model. The point interactions will 
be treated in the standard way || . We suppose that they 
are finitely many and supported by a set Y := { yj : j — 
1, . . . , N } . In the vicinity of the point yj any solution 
to the stationary Schrodinger equation behaves as 
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where the coefficients are related by Bj 
and the real parameter ctj characterizes the interaction 
"strength" §. 

Suppose that the incident particle momentum is par- 
allel to the first axis and equals k . For a given family of 
the coupling constants, a :— {ct\, . . . , a^} , the scatter- 
ing wavefunction equals [pi Sec. II. 1.5] 
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with x = (x, y, z) , yj = {y jl ,y j2 ,y J 3) , and 



T a , Y (k) 
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where Gk is the regularized Green's function, 

Gk(x) = - 

47r|x| 

for \x\ ^ and zero otherwise. 

We denote by 7Z, T the real and imaginary part of the 
function f/' Q ,y(fc; •) , respectively; to find nodal sets one 
has to solve the equations 



fc(x,y,z) = = l(x,y,z) . 



(4) 



Some conclusions about the existence and properties of 
the solutions can be derived from the implicit-function 
theorem 0: 

(i) In the absence of the scatterers, solutions to each 
of the conditions (0) are planes perpendicular to 
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the x-axis. They can be used as "unperturbed" 
solutions, since 

S(( n+ £) 

as y 2 + z 2 — ■> oo . Hence far enough from the scat- 
terers the solutions to the two conditions are locally 
unique. 

(ii) Whenever a unique solution exists, it is locally a 
C°° surface, since TZ, X are real-analytic functions 
of the variables x,y,z. For the same reason, the 
sought nodal sets are C°° curves being intersec- 
tions of two surfaces, except for possible crossing 
points where the solution is locally nonunique. 

(Hi) The implicit-function theorem also implies that the 
nodal lines are confined to a bounded region of the 
configuration space only. Indeed, the zero surfaces 
are for large y 2 + z 2 of the form 

xi = (n+ n + f n (y,z), x x = nir + g n (y , z) , 

respectively, and the smooth functions /„, g n are 
O ((y 2 + z 2 )" 1 / 2 ) , so that there is no intersection 
for y 2 +z 2 large enough; for large \n\ the argument 
can be extended up to the x-axis. 

Before we shall discuss whether the solutions do exist, let 
us point out two other properties: 

(iv) It follows from (Hi) that the nodal line, if they exist, 
consist of a family of closed loops, each of them 
being a C°° curve. They can cross; this happens if 
one of the surfaces has locally a saddle shape and 
the other one — apropriately shifted by parameter 
choice — is "flatter", i.e., has a small enough mean 
curvature. 

(v) On the other hand, the loops cannot entagle into 

nontrivial knots, because any closed loop on a 
smooth surface without intersections is topologi- 
cally equivalent to a circle; one cannot "draw" a 
knot upon a surface. 

Before proceeding further, let us remark that the above 
argument has two basic ingredients. One is the finite 
range of the interaction which implies that far from the 
scattering centers the zero surfaces of the functions TZ 
and T cannot intersect. The other is the analyticity of 
the wavefunction with respect to the coordinates. Since 
this property persists if point scatterers are replaced by 
finite hard obstacles or a smooth enough potential, the 
conclusions are expected to apply to a much wider class 
of three-dimensional scattering systems. 

To show that the described scheme is not empty, let 
us discuss next the case of a single 6 potential of the 
"strength" a situated at the origin. In view of the cylin- 
drical symmetry, possible nodal lines are circular loops 
perpendicular to the x-axis, hence we replace (y, z) by 



(y cos tp^y sin <p) so |a?| = x 2 + y 2 . The scattering so- 
lution (|2|) now acquires the form 



^ a fi{k-x) = e lkx + 



ik\x\ 



(Ana — ik)\x\ 
where the right side turns to zero if 

e ik(\x\-x) ik 
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Taking the modulus, we find the distance of the ring from 
the origin, 
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y/k 2 + 167T^Q 



2„,2 



(7) 



Furthermore, comparing the real and imaginary part in 
(H) we find 



1 a * 
— Arctan 

k Ana 
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where the arctangent branch remains to specified. De- 
noting the right side of (H) as —7 and substituting into 



x 2 + y 2 , we get 
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7 2 — 27a: = —7(7 + 2|x|) . 



(9) 



A nontrivial nodal ring exists if the left side is positive. 
It is straightforward to see from here that this happens if 
— 2\x\ < 7 < 0; introducing k :— k/Aira we can rewrite 
the last condition as 



2|/s| 



< Arctan k < . 



(10) 



At the edges of the interval the ring shrinks into a nodal 
point. Inspecting the inequalities (|l0|), we infer that for 
a < there is a unique solution given by the basic 
branch. On the other hand, one has to use arctan k — wk 
if a > and the solution exists provided k > 2.971 , i.e. 



- < 2.679 x 10" 

k 



(11) 



If there is more than one point scatterer, it is no longer 
possible to find the nodal lines analytically. However, 
the explicit form ^ of the scattering wavefunction al- 
lows a numerical treatment. This demonstrates not only 
that the nodal loops can exist in the multicenter situation 
too, but also their possible crossing described in point (v) 
above; changing appropriately the parameters, one can 
achieve that a loop undergoes a fission process in which it 
is first "strangled" and then decays into two loops. As an 
illustration we show in Fig. 1. nodal lines inside a given 
rectangular "box" corresponding to 10 scattering centers 
whose positions are marked by dots. All the point inter- 
actions have the same strength a = . To convey the 
their spatial shapes, we do not plot the lines themselves 
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but rather their properly lighted tubular neighborhoods 
of radius R = 0.1 . 

To understand the behavior of the wavefunction in the 
vicinity of the nodal lines, one has to realize that due to 
the smoothness of the latter, the wavefunction has locally 
an approximate cylindrical symmetry. Hence if we want 
to expand i]) a ,Y(k\ x) , it is natural to place the point of 
interest into the origin with the nodal line tangent to the 
z axis, we can use again the coordinates of the above 
example and write 

*p a ,Y(k;r) = Y, c ^ m<p J m (ky). (12) 

The m = term is missing because the wavefunction 
vanishes by assumption at y = . The phase behavior 
around the node is determined by the terms with the low- 
est nonzero \m\ . Suppose that c\ or c_i is nonzero and 

|ci| 7^ | c x | • Up to higher-order terms, the probability 

current j(f) is then located in a plane perpendicular to 
the nodal line and 

J^v{r)dr = 27rsgn(|ci| - |c_i|) ; (13) 

where v denotes the velocity of the probability flow, 
v(r) :— V<f>(r) ■ By analyticity this result extends to 
any curve C encircling the nodal line once as long as 
it stays away of the scattering centers. It is moreover 
clear that the decribed case is generic, because it cor- 
responds to simple zeros of the wavefunction. Higher 
winding numbers may appear only if there is too > 
such that c m = for m = — m , — m + 1, . . . , to . 

To illustrate these considerations we choose a part of 
the "lowest" nodal loop of Fig. 1. To show that this is 
indeed the |m| = 1 case, we plot in Fig. 2 the wave- 
function phase modulo 2tt in the planes perpendicular 
to the loop. A single cut is clearly visible. In Fig. 3 we 
plot currents vectors at a tubular surface centered at the 
nodal line. 

Let us ask whether the existence of "tornados" can 
have observable consequences. We have remarked that 
scattering of charged particles may produce in this way 
closed magnetic flux lines. Consider therefore a scatter- 
ing of electrons and heavy ions in the presence of fixed 
obstacles. If the ion-obstacle configuration is properly 
chosen, the electrons should exhibit interference behav- 
ior of the AB type. 

In conclusion, we have shown here analytically and 
demonstrated numerically that the probability flux in 
three-dimensional scattering can exhibit vortical behav- 
ior around nodal lines of the wavefunction. The latter 
are generically closed smooth curves which may intersect 
but do not entangle into knots, with the probability cur- 
rent velocity integral over a curve encircling the nodal 
line once equal to ±27r . 



ACKNOWLEDGMENTS 

The research has been partially supported by the The- 
oretical Physics Foundation in Slemeno and by GACR 
under the contract 202-0218/96-98. 



[1] P.A.M. Dirac, Proc. Roy. Soc. A133, 60 (1931); A. A. 
Abrikosov, JETP 32, 123 (1957); Y. Aharonov, D. Bohm, 
Phys. Rev. 56, 485 (1959). 

[2] For a rigorous analysis of GL vortices and further refer- 
ences see F. Bethuel, H. Brezis, and F. Helein, Ginzburg- 
Landau Vortices, Birkhauser 1994; E.H. Lieb and M. Loss, 
Math. Res. Lett. 1, 701 (1994); Yu.M. Ovchinnikov, 
I.M. Sigal, Ginzburg-Landau Equation, Preprint Toronto 
1997. 

[3] K.-F. Berggren, C. Besev, Zhen-Li Ji, Physica Scripta 
T42, 141 (1992); K.-F. Berggren, Zheng-Li Ji, Phys. 
Rev. B47, 6390 (1993); S. Chaudhuri, S. Bandyopad- 
hyay, M. Cahay, Phys. Rev. B45, 11126 (1992); P Exner, 
P. Seba, M. Tater, D. Vanek, J. Math. Phys. 37, 4867 
(1996); V.M. Ramaglia, F. Ventriglia, GP Zucchelli, 
Phys. Rev. B48, 2445 (1993); B52, 8372 (1993); Hua 
Wu, D.W.L. Sprung, Phys. Lett. A183, 413 (1993); 
K.N.Pichugin, A.F.Sadreev, JETP 82, 290 (1996). 

[4] Since values of the constants are not important in the 
present letter, we put h = 2m = 1 . 

[5] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, H. Holden, 
Solvable Model in Quantum Mechanics, Springer 1988. 

[6] The free case corresponds in this formalism to a — oo . 
Approximation of these S interactions in terms of scaled 
potentials is discussed in detail in Ref. jEj. If ay < one 
can interpret such a potential also as a low-energy de- 
scription of a narrow nonresonant potential well or barrier 
of the radius (— A-KUj)" 1 — cf. P. Exner, P. Seba, Phys. 
Lett. A222, 1 (1996). 

[7] L. Schwartz, Analyse Mathematique, Hermann: Paris 
1967, Sec. III. 8; V. Jarmk, Differential Calculus II, 
Academia: Prague 1956, Sec. VIII. 1. 



Figure captions 

Figure 1. Nodal lines of the wavefunction (shown by 
means of their lighted tubular neighborhoods) for 
scattering on 10 6 potentials of the same strength 
a = , whose positions are marked by dots, with 
the incident particle momentum k = 2 . 

Figure 2. Wavefunction phase around the "lowest" 
loop of the preceding picture. The wavefunction 
phase </>(F) has been evaluated mod 27r on a tubu- 
lar surface of radius 0.1 along the nodal line and 
plotted in the perpendicular planes being scaled 
down by the factor 0.02. 
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Figure 3. The probability flow at the same tubular sur- 
face as in Fig. 2, the observation point is changed. 
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